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We discuss the relaxation kinetics of a one-dimensional dimer adsorption model as recently 
proposed for the binding of biological dimers like kinesin on microtubules. The non-equilibrium 
dynamics shows several regimes: irreversible adsorption on short time scales, an intermediate plateau 
followed by a power-law regime and finally exponential relaxation towards equilibrium. In all four 
regimes we give analytical solutions. The algebraic decay and the scaling behaviour can be explained 
by mapping onto a simple reaction-diffusion model. We show that there are several possibilities to 
define the autocorrelation function and that they all asymptotically show exponential decay, however 
with different time constants. Our findings remain valid if there is an attractive interaction between 
bound dimers. 



Motor proteins are of fundamental importance for intracellular transport and many other biologically relevant 
transport processes. Recently, a large body of data on these proteins has been collected using a diverse set of 
experimental tools ranging from single-molecule mechanics to biochemical methods Q|. In many instances these 
experimental systems are both versatile experimental techniques and interesting non-equilibrium model systems. One 
system of particular interest is a standard method from biophysical chemistry known as "decoration experiments" . 
Here monomeric or dimeric motor enzymes are deposited on their corresponding molecular tracks at high densities (see 
fig. |l|a). For kinesin motors, these tracks are microtubules, hollow cylinders usually consisting of 13 protofilaments. 
The kinesin binding sites are located on the P subunits which form a helical (wound-up rhombic) lattice with a 
longitudinal periodicity of 8 nm. Decoration techniques have traditionally been used to investigate the structure 
and the binding properties of kinesin i. e. after waiting for the system to equilibrate the binding stoichiomctry is 
determined and the structure by cryo-electron microscopy followed by 3D image reconstruction. 

These experiments call for a theoretical analysis of dimer adsorption kinetics with competing single and double 
bound dimers and a finite detachment rate. For a quantitative analysis of decoration data ^ one needs to know 
the binding stoichiomctry in the equilibrium state in terms of binding constants for the first and second head of the 
dimer molecule. The dynamics of the approach to equilibrium is useful to estimate when an experimental system can 
be considered as equilibrated Q. Even more importantly, time- resolved decoration experiments (e.g. by using motor 
enzymes labelled by some fluorescent marker) combined with our theoretical analysis could provide new information 
about reaction rates which are to date not known completely. Understanding the kinetics of passive motors is 
undoubtedly a necessary prerequisite for studying the more complicated case of active motors at high densities. The 
model is also interesting in its own right since it contains ~ as detailed below - some novel features of non-equilibrium 
dynamics of dimer models which have not been discussed previously. 

There seems to be convincing evidence that kinesin heads can bind on two adjacent binding sites only in longitudinal 
but not in lateral direction ||^ . This introduces a strong uniaxial anisotropy and distinguishes the adsorption process 
of protein dimers from simple inorganic dimers. If we take into account only steric interactions and (for now) neglect 
nearest neighbour attractive interaction, we are left with a one-dimensional problem of kinesin dimers decorating 
a single protofilament (one-dimensional lattice). Then our model is defined as follows. Kinesin is considered as a 
dimeric structure with its two heads tethered together by some flexible joint. Hence each dimer (kinesin protein) can 
bind one of its two heads (motor domains) to an empty lattice site The binding rate fc+i c for this process is 
proportional to the solution concentration c of the dimeric proteins. Successively, the dimer may either dissociate from 
the protofilament with a rate fc_i or also bind its second head to an unoccupied site in front (f) of or behind (b) the 
already bound head. Since kinesin heads and microtubules arc both asymmetric structures the corresponding binding 
rates ^-^^ ^+2 general different from each other. The reverse process of detaching a front or rear head 

occurs at rates A:;^2 ^^'^ ^-2- ^ reaction scheme with all possible processes and their corresponding rate constants is 
shown in fig. |^b. 



1 




Fig. §b. Fig. I 

FIG. 1. a) Kinesin dimers can bind one head or both heads to a one-dimensional lattice (tubulin protofilament). Binding 
sites are located on /3-tubulin subunits (dark) while a-subunits (bright) are irrelevant for our study, b) Reaction scheme for all 
binding and unbinding processes. 

FIG. 2. Vacancy concentration as a function of time t for K = k+c/k- = 1, ... , 10000. Dashed lines show the steady state 
concentration no = 1/ \/l + at given K. The line with circles shows the short-time limit for K = 10000 |5| , and the dashed 
line the result of the reaction-diffusion model (Eq. ^). The thin dot-dashed curve shows the behaviour for K = 100 with an 
attractive interaction between bound dimers with parameters yl = 10 and B = 1/10 (a dimer is 10 times as likely to associate 
to and 10 times less likely to dissociate from a certain site if one of the neighbours is present) . 

Since there is no energy source in the system (unlike situations where ATP is being hydrolysed), the reaction rates 
obey the principle of detailed balance. As a consequence, the reaction rates are actually not independent from each 
other but the ratio of the on- and off-rates has to equal the equilibrium binding constant K2 — = 
similarly we have Ki = k^i/k^i. A particular coverage of the lattice is described as a sequence of dimers bound 
with both heads (D), one head only (1) and empty sites (0). We denote the probabilities to find a certain lattice 
site in one of these states by 2n£), ni and ng, respectively. Of course, normalisation of the probabilities requires 
no + ni + 2n_D = 1. 

We start our analysis with the stoichiometry of the final steady state. Important information can be gained already 
by exploiting the symmetries of the kinetic process. The reaction scheme in fig. ^3 (including steric constraints) 
remains invariant under the following transformation: (I) all dimers bound on a single head are replaced by vacant 
sites and vice versa (0 ^ 1), while dimers bound with both heads are kept as they are {D ^ D). (II) The attachment 
and detachment rate of the first head are exchanged and the forward/backward binding rates of the second head are 
mirrored (fc-i^ic, fc^j) ^ (^-ii fc±2)- Since the coverage in the steady state is only a function of Kic and K2 (see below) 
this symmetry implies that the mean number of dimers attached with both heads no{Kic,K2) is invariant against 
interchanging the attachment and detachment rates of the first head, njj^l/ (Kic), K2) — nD{Kic, K2). Similarly, the 
mean total number of bound heads per lattice site (often simply called binding stoichiometry), v ~ 2{ni +n£i), obeys 
the symmetry relation i>{l/{Kic), K2) = 2 — v{Kic, K2). Thus we conclude that reaches its maximum at Kic = 1 
where v = 1. 

We determine the actual value of the mean occupation numbers in the steady state using detailed balance and 
the fact that the dimers have only a hard-core interaction. Hence the probability to find a certain sequence of 
O's, I's and D's {e.g. "0,1,D,D,1,D") has to be invariant against permutations of these states. In such a random 
sequence the probabilities to find a particular state 0, 1 or _D at a certain place are given by pi — rii/ {no + ni + no). 
Detailed balance requires that for each pair of possible configurations, their probabilities are in the same ratio as 
the transition rates between them. Hence the ratio of probabilities to find a sequence with a 1 or at a certain 
place is Pi/pq — fc+ic/fc-i = Kic. Similarly, we get for transitions between D and 01: po/poPi = ^+2/^-2 ~ -^'2- 
These two equations, together with the normalisation condition, uniquely determine the values no, ni and ud. The 
stoichiometry, i.e. the total number of heads per binding site = 2(ni + nn) is given by 



2 



The number of dimers bound with both heads per lattice site reaches its maximum = (1 — Xj^/K^ + l)/2 for 

We now turn to the dynamics of the biological dimcr model. In order to avoid unnecessary complications we restrict 
ourselves to the limit K2 — > 00 with K = K1K2C fixed. Then our model reduces to a dimer deposition-evaporation 
model where the dimers can only bind and unbind with both heads at the same time. The vacancy concentration 
in the steady state then simplifies to rig — 1/ Vl + '^K. This simplified model still captures all the essential aspects 
of the non-equilibrium dynamics. Similar dimer adsorption models have been studied previously |^^. Privman and 
Nielaba [|| studied the effect of diffusion on the dimer deposition process. There are several key differences to our 
model, the most important of which is that diffusion without detachment results in a 100% saturation coverage, 
whereas a model with detachment leads to a limiting coverage whose value depends on the binding constants of the 
first and second head (see eq. |l|). This also has important implication on the dynamics as discussed below. For 
example, as a consequence of a finite coverage the final approach to equilibrium is not a power law but exponential 
and there are, in addition, interesting temporal correlations in the fluctuations in the steady state. Stinchcombe and 
coworkers 1^,^ studied the effect of detachment on the adsorption kinetics but allowed for regrouping of attached 
dimer molecules (two monomers that belonged to different dimers during attachment can form a dimer and detach 
together). While such processes are allowed for some types of inorganic dimers, they are certainly forbidden for dimer 
proteins like kinesin where the linkage between its two heads is virtually unbreakable. If regrouping is allowed the 
steady state auto-correlation functions for the dimer density shows an interesting power-law decay cx If if 

is forbidden this power-law decay is lost (and becomes an exponential to leading order) due to the permanent linkage 
between the two heads of the dimer. Intuitively this may be understood as follows: only if regrouping of dimers is 
allowed are there locally jammed configurations (Neel-like states, "101010" , with alternating occupied and unoccupied 
sites in which neither attachment nor detachment of dimers is possible) in the final steady state which slow down 
the dynamics. In addition, we will show that the autocorrelation functions of the dimer and the vacancy occupation 
number show strong differences in shape and typical times scales of relaxation. 

The basic kinetic steps in the reduced model are deposition without overlap and evaporation without regrouping 
of the dimers; the effective attachment and detachment rates k± in terms of the original model are k± — ki±{k2± + 
k\^)/{k^i -t- ki^2 + ^+2)- Processes in which one head detaches on one side and subsequently attaches on the other 

side also lead to explicit diffusion with a rate — fc^2^+2/ (^2(^-1 + ^+2 "I" ^+2))- This feature is not essential since 
it can be incorporated in the effective diffusion introduced later. It will hence be disregarded in what follows (if one of 
the rates fc|.2 is small, it is negligible anyway). Note that as a consequence of detailed balance the diffusion of dimers 
is symmetric despite the asymmetry in reaction rates. 

To study the kinetics we choose the initial condition as typically used in an experiment, namely an empty lattice. 
Figure ^ shows simulation data Q for the average vacancy concentration as a function of time for a set of binding 
constants K = k^c/k^ = K1K2C. We find qualitatively very different approaches to the final steady state depending 
on the value of the constant K. For if <C 1, where the off-rate fc_ is much larger than the on-rate fc+c, there is 
no crowding on the lattice and the dimeric nature of the molecules does not affect the approach to equilibrium, 
which is, like for monomers, exponential with a decay rate In the opposite limit, if ^ 1, we find a two- 

stage relaxation towards the steady state. The vacancy concentration as a function of time reveals four regimes, an 
initial attachment phase, followed by an intermediate plateau, then a power-law decay and finally an exponential 
approach towards equilibrium. At short time scales, t <C kZ^ , when only deposition processes are frequent but 
detachment processes are still very unlikely, the kinetics of the model is equivalent to Flory's random sequential 
dimer adsorption model [ pO| , pT| . There is an initial decay obeying the kinetics described in |^,|l^ with a vacancy 
concentration no(t) — exp(— 2 + 26^*^+^*) . The vacancy concentration locks at an intermediate plateau ng = 
p0{ in the time interval between the characteristic attachment time t_|_ = l/(fc_|_c) and detachment time r_ = l/fe_. 
This Flory plateau represents a configuration in which all remaining vacancies are isolated, causing the system to be 
unable to accommodate for the deposition of additional dimers. The secondary relaxation process towards the final 
steady state is enormously slowed down. It shows a broad time domain with a power-law oc t~^/'^ instead of a simple 
exponential decay. Similar multi-stage relaxation processes have been observed in dimer models with diffusion but no 
detachment |^ , the key difference being that the detachment process implies that the steady state has a finite vacancy 
density and the final approach to the steady state remains not a power law but becomes exponential. There are also 
interesting similarities. In particular, in both models a large portion of the final approach to the steady state is 
mediated by the annihilation of vacancies. This behaviour can be explained by introducing a particle representation 
in the following way (analog to the adsorption-diffusion model |^). We denote each vacancy on the lattice as a 
"particle" A, and each bound dimer as an inert state (00). 
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FIG. 3. Time evolution of a state, consisting of attachment (fast) and detachment (slow) events. The coarse-grained inter- 
pretation includes only long-living states (in boxes). The effective steps include diffusion, pair annihilation and pair creation 
(not shown). 



FIG. 4. Scaled equilibrium autocorrelation functions for Co{t) = KCnit), Cn(t) — Coit) and Cp{t) — KCp{t) , t — k-t/K. 
The function Cp{i) is compared with the analytical calculation from ref. |0 (solid line). Simulation data were obtained with 
K = 100 (circles, for Co and Cd connected with dot-dashed and dashed lines) and K = 400 (triangles) . 



The detachment of a dimer then corresponds to a pair creation process 00 — *■ AA, and the decoration process to 
pair annihilation AA — > 00. Since we consider the limit K ^ 1, states with two vacancies A on neighbouring sites 
have a very short lifetime. We may thus introduce a coarse-grained model by eliminating these states (see fig. ||). 
Then processes like AOO AAA 00 A result in an effective diffusion for particle A with a hopping rate rhop = fc_/2 
and an effective step width of two lattice sites. Pair annihilation, AA 00, occurs with a rate fc+c. Pair creation, 
on the other hand, occurs mainly through the process 0000 OOAA AAAA AOOA; the corresponding rate is 
k-/3K per lattice site and hence largely suppressed with respect to the annihilation process as long as the particle 
concentration is far above its steady-state value. Other processes involve more particles and are of higher order in 
terms of a power series in K~^. They are negligible for ua <C 1 and K ^ 1. In summary, for K ^ 1, our dimer 
model can be mapped onto a one-particle reaction-diffusion model A + A ^ 0. Pair creation processes, ^ A + A, 
are highly suppressed and do not play any role until the system comes close the steady state. 

Models like the simple reaction-diffusion model A + A —> show interesting non-equilibrium dynamics [p^|-p^ . 
One can show ||l^,|l8| that asymptotically the particle density n^(i) decays algebraically, n^(t) oc t^^/^ which nicely 
explains the slow decay observed in simulation data (see fig. |^). Note that a mean-field like rate equation approach 
would predict riAit) cx t~^. In our analysis we can even go beyond the asymptotic scaling analysis and try to compare 
with exact solutions of the model for a random initial distribution with density p by Krebs et al. [ p^ . They find 
(adapted to our situation with two-site hopping) 



du 



(2) 



Its asymptotic limit (first determined by Torney and McConnel ||T^) reads ua = (327rrhopi)^^^^; note that it is 
independent of the initial particle concentration in the Flory plateau p. 

Our Monte-Carlo data (see fig. ^ are in excellent agreement with the predictions of eq. |^. Minor deviations at 
times between the plateau and the power-law decay are due to the assumption of a random particle distribution 



underlying the derivation of ec 
intermediate plateau regime |1^ 



3; the random sequential adsorption process leads to some particle correlations in the 
[17 but they do not affect the asymptotic behaviour as they are only short-ranged. The 
results from the ^ -I- ^ ^ model also become invalid for very long times where the particle concentration comes 
close to its equilibrium value. In this limit the dynamics becomes scale-invariant pQ]. The particle concentration can 



4 



be written in scaling form nA{t) = K~^^'^n{r-iiopt / K) with a diverging characteristic time scale tk oc K. 

A limitation of the above mapping becomes evident if one considers the equilibrium autocorrelation functions. 
Contrary to conventional models there are three different autocorrelation functions with different functional forms 
(fig. ^). Co(t) = {fiOi{to)fi'Oi{to + 1)} — ("-0;)^ describes the correlation function of the probability to find a vacancy 
at a certain lattice site. C^i is the equivalent quantity defined from the probability to find a dimer on a certain 
pair of sites hjj. .^^ and Cp from the probability to find a vacancy on at least one from a pair of neighbouring 
sites, ho- + ?T-Oi+i — no-ho.^-^. For K ^ I the autocorrelation functions become scale invariant as well. Their scaling 
form reads Co{i) = KCo{t), Coii) = Coit) and Cp{i) = KCp{t) with i = k_t/K. The latter corresponds to the 
autocorrelation function in a reaction-diffusion model A + A 0, which has recently been calculated analytically in 
ref. (l| 

Cp{t) = -= - Erfc V 2t Erfc \/2t , with Erfc (x) = -= / dy . (3) 

\ \/9.7rf / VTT . L 



The other two functions decay on the same time-scale, but with different prefactors. The reason is that even if a pair 
of vacancies annihilates, the system still keeps memory on whether the surrounding dimers were located on even or 
odd locations and this gives those correlation functions that distinguish between even and sites a longer decay time. 

Finally, we would like to note that an attractive interaction between attached dimers plays an important role in 
some cases. In the case of kinesin, there has been an observation of coexisting empty and decorated domains which can 
only be explained by an attractive interaction Similar observation has been done on actin decorated with myosin 
pl| and tropomyosin . We introduce the interaction by assuming that a dimer is more likely to bind to a pair sites 
if one or two neighbours are already bound. The binding rate then becomes Ak+ (one neighbour bound) or 
(both neighbours bound). Similar, we assume that a dimer with one bound neighbour dissociates with rate _Bfc_ and 
that a dimer with two bound neighbours dissociates with rate B^k^. This interaction changes both relaxation stages 
quantitatively. First, the vacancy concentration on the intermediate plateau lowers since the interaction improves 
the formation of contiguous clusters during the first stage. Second, the diffusional relaxation slows down since the 
detachment rate decreases. And finally, the equilibrium vacancy concentration decreases. Nevertheless, interacting 
models show the same two-stage relaxation behaviour. An example of a model with interaction is shown by the 
dot-dashed line in fig. ^ We therefore expect that many qualitative conclusions from this article will apply to a much 
wider range of biological adsorption problems such as the binding of double-headed myosin or tropomyosin . 
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